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Abstract 

A search for excited neutrinos is performed using the full e~p data sample collected by 
the HI experiment at HERA at a centre-of-mass energy of 319 GeV, corresponding to a to- 
tal luminosity of 184 pb^^. The electroweak decays of excited neutrinos i/*^i^7, i'*^i'Z 
and u*^eW with subsequent hadronic or leptonic decays of the W and Z bosons are con- 
sidered. No evidence for excited neutrino production is found. Mass dependent exclusion 
limits on u* production cross sections and on the ratio of the coupling to the compositeness 
scale //A are derived within gauge mediated models. A limit on //A, independent of the 
relative couplings to the SU(2) and U(l) gauge bosons, is also determined. These limits 
extend the excluded region to higher masses than has been possible in previous excited 
neutrino searches. 
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1 Introduction 



The three-family structure and mass hierarchy of the known fermions is one of the most puz- 
zling characteristics of the Standard Model (SM) of particle physics. Attractive explanations 
are provided by models assuming composite quarks and leptons [1]. The existence of excited 
states of leptons and quarks is a natural consequence of these models and their discovery would 
provide convincing evidence of a new scale of matter. Electron-proton interactions at very high 
energies provide a good environment in which to search for excited states of first generation 
fermions. In particular, excited neutrinos (u*) could be singly produced through the exchange 
of a boson in the t-channel. 

In this letter a search for excited neutrinos using the complete e~p HERA collider data of 
the HI experiment is presented. Electroweak decays into a SM lepton (e, z/g) and a SM gauge 
boson (7, W and Z) are considered and hadronic as well as leptonic decays of the W and Z are 
analysed. 

The data, collected at electron and proton beam energies of 27.6 GeV and 920 GeV, respec- 
tively, correspond to a total integrated luminosity of 184 pb"^. During the HERA II running 
period, the electron beam was longitudinally polarised, at a level of typically 35%. For this 
analysis the periods with left-handed and right-handed beams are combined and the analysed 
data sample has a residual polarisation of 5% left-handed. With more than a ten-fold increase 
in statistics, this analysis supersedes the result of the previous HI search for excited neutrinos 
based on a data sample corresponding to a luminosity of 15 pb^ [2]. 



2 Phenomenology 



In the present study a model [3-5] is considered in which excited fermions are assumed to 
have spin 1/2 and isospin 1/2. Both left-handed and right-handed components of the excited 
fermions form weak iso-doublets F£ and F^. In order to prevent the light leptons from radia- 
tively acquiring a large anomalous magnetic moment [6,7], only the right-handed component 
of the excited fermions takes part in the generalised magnetic de-excitation. The interaction be- 
tween excited fermions, gauge bosons and ordinary fermions is then described by the effective 
Lagrangian [4]: 



int. 



2A ^ 



^11/ 



A" 



Fl + h.c. 



(1) 



The matrix a^" is the covariant bilinear tensor, VF^^, B^i, and G^^, are the field-strength 
tensors of the SU(2), U(l) and SU(3)c gauge fields, r", Y and A" are the Pauli matrices, the 
weak hypercharge operator and the Gell-Mann matrices, respectively. The standard electroweak 
and strong gauge couplings are denoted by g, g' and g^, respectively. The parameter A has 
units of energy and can be regarded as the compositeness scale which reflects the range of a 
new confinement force. The constants /, /' and fg are form factors associated to the three 
gauge groups. They can be interpreted as parameters setting different scales Aj = A//j for the 
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different gauge groups, thus allowing the composite fermion to have arbitrary coupling strengths 
with the three gauge bosons. 

Following this model, single production of excited neutrinos in ep collisions may result from 
the t-channel exchange of a FT boson. Due to the helicity dependence of the weak interaction 
and given the valence quark composition and density distribution of the proton, the v* produc- 
tion cross section is predicted to be much larger for e^p collisions than for e^p. For a u* mass 
My. of 200 GeV the ratio of the production cross sections is of order 100. The u* production 
cross section is expected to scale linearly with the polarisation of the incident electron beam, 
similarly to the SM charged current process. The excited neutrino may decay into a lepton 
and an electroweak gauge boson via u* — >z/7, u*^eW and ij*~^ijZ. For a given My* value 
and assuming a numerical relation between / and /', the v* branching ratios are fixed and the 
production cross section depends only on //A. The u* is expected not to have strong interac- 
tions and therefore this search is insensitive to fg. Two complementary coupling assignments 
/ = +/' and / = — /' are studied in detail. For / = +/' the excited neutrino has no tree-level 
electromagnetic coupling [8] and therefore the photonic decay of the v* is forbidden whereas 
for / = — /' decays into z/7, uZ and eW are allowed. In addition, arbitrary ratios of /'// are 
considered in the range —5 to +5. 



3 Simulation of Signal and Background Processes 

A Monte Carlo (MC) program developed for this analysis is used for the calculation of the u* 
production cross section and the simulation of signal events. The events are simulated using 
the cross section calculated from the Lagrangian described in equation ([T]) using the CompHEP 
package [9]. Initial state radiation from the incident electron is included using the Weizsacker- 
Williams approximation [10]. The proton parton densities are taken from the CTEQ5L [11] 
parametrisation and are evaluated at the scale \/Q^, where Q"^ is the four-momentum transfer 
squared. The parton shower approach [12] is applied to simulate Quantum Chromodynamics 
(QCD) corrections in the initial and final states. The hadronisation is performed using Lund 
string fragmentation as implemented in PYTHIA [12]. 

In the MC generator the full transition matrix including the production and the decay is 
implemented. This is important if the natural width of the u* is large, which is typically the case 
at high mass where factorisation of the u* production and its decay no longer holds. Events used 
in the determination of signal efficiencies are generated with the coupling //A corresponding, 
for each u* mass, to the expected boundary of the probed domain in the plane defined by M^* 
and //A. 

The Standard Model background processes that may mimic a u* signal are neutral current 
(NC) and charged current (CC) deep-inelastic scattering (DIS) and to a lesser extent photopro- 
duction, lepton pair production and real W boson production. 

The RAPGAP [13] event generator, which implements the Bom, QCD Compton and Boson 
Gluon Fusion matrix elements, is used to model NC DIS events. The QED radiative effects 
arising from real photon emission from both the incoming and outgoing electrons are simulated 
using the HERACLES [14] program. Direct and resolved photoproduction of jets and prompt 
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photon production are simulated using the PYTHIA event generator. The simulation is based 

on Bom level hard scattering matrix elements with radiative QED corrections. In RAPGAP and 
PYTHIA, jet production from higher order QCD radiation is simulated using leading logarith- 
mic parton showers and hadronisation is modelled with Lund string fragmentation. The leading 
order MC prediction of NC DIS and photoproduction processes with two or more high trans- 
verse momentum jets is scaled by a factor of 1.2 in order to normalise to next-to-leading order 
QCD calculations [15]. Charged current DIS events are simulated using the DJANGO [16] pro- 
gram, which includes first order leptonic QED radiative corrections based on HERACLES. The 
production of two or more jets in DJANGO is accounted for using the colour-dipole-model [17]. 

Contributions from elastic and quasi-elastic QED Compton scattering are simulated with 
the WABGEN [18] generator. Contributions arising from the production of W bosons and 
multi-lepton events are modelled using the EPVEC [19] and GRAPE [20] event generators, 
respectively. 

All processes are generated with an integrated luminosity significantly higher than that of 
the data sample. Generated events are passed through the full GEANT [21] based simulation of 
the HI apparatus, which takes into account the running conditions of the different data taking 
periods, and are reconstructed and analysed using the same program chain as for the data. 



4 Experimental Conditions 

A detailed description of the HI experiment can be found in [22]. Only the detector compo- 
nents relevant to the present analysis are briefly described here. The origin of the HI coordinate 
system is the nominal ep interaction point, with the direction of the proton beam defining the 
positive 2;-axis (forward region). Transverse momentum (Pt) is measured in the xy plane. The 
pseudorapidity r] is related to the polar angle 6^ by 77 = — In tan(^^/2). The Liquid Argon (LAr) 
calorimeter [23] is used to measure electrons, photons and hadrons. It covers the polar angle 
range 4° < ^ < 154° with full azimuthal acceptance. Electromagnetic shower energies are 
measured with a precision of a{E)/E = 12%/^£;/GeV © 1% and hadronic energies with 
a(E)/E = 50%/^/E/GeV © 2%, as measured in test beams [24,25]. In the backward re- 
gion, energy measurements are provided by a lead/scintillating-fiber (SpaCal) calorimeter [26] 
covering the angular range 155° < < 178°. The central (20° < 9 < 160°) and forward 
(7° < 6* < 25°) tracking detectors are used to measure charged particle trajectories, to re- 
construct the interaction vertex and to complement the measurement of hadronic energy. The 
LAr and inner tracking detectors are enclosed in a super-conducting magnetic coil with a field 
strength of 1.16 T. The return yoke of the coil is the outermost part of the detector and is 
equipped with streamer tubes forming the central muon detector (4° < < 171°). In the for- 
ward region of the detector (3° < 6* < 17°) a set of drift chambers detects muons and measures 
their momenta using an iron toroidal magnet. The luminosity is determined from the rate of 
the Bethe-Heitler process ep— >ep7, measured using a photon detector located close to the beam 
pipe at z — —103 m, in the backward direction. 
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5 Data Analysis 



The triggers used in this analysis are based on the detection of energy deposits in the LAr 
calorimeter [27]. Events containing an electromagnetic deposit (electron or photon) with an 
energy greater than 10 GeV are triggered with an efficiency close to 100%. For events with 
missing transverse energy above 20 GeV, the trigger efficiency is ~ 90%. 

In order to remove background events induced by cosmic showers and other non-ep sources, 
the event vertex is required to be reconstructed within 35 cm in z of the nominal interaction 
point. In addition, topological filters and timing vetoes are applied. 

The identification of electrons or photons relies on the measurement of a compact and iso- 
lated electromagnetic shower in the LAr calorimeter. In addition, the hadronic energy within 
a distance in the pseudorapidity-azimuth (77 — 0) plane R — y^A^^~+A0^ < 0.5 around the 
electron (photon) is required to be below 3% of the electron (photon) energy. Muon iden- 
tification is based on a track measured in the inner tracking systems associated with signals 
in the muon detectors [28]. A muon candidate should have no more than 5 GeV deposited 
in a cylinder, centred on the muon track direction, of radius 25 cm and 50 cm in the elec- 
tromagnetic and hadronic sections of the LAr calorimeter, respectively. Calorimeter energy 
deposits and tracks not previously identified as electron, photon or muon candidates are used to 
form combined cluster-track objects, from which the hadronic energy is reconstructed [29,30]. 
Jets are reconstructed from these combined cluster-track objects using an inclusive kx algo- 
rithm [31, 32] with a minimum transverse momentum of 2.5 GeV. The missing transverse mo- 
mentum P™^^ of the event is derived from all identified particles and energy deposits in the 
event. The P™'**** is assumed to originate from a single neutrino. The four-vector of this neu- 
trino candidate is reconstructed assuming transverse momentum conservation and the relation 
Y.i{E'' - P^) + {E" - P^) = 2El = 55.2 GeV, where the sum runs over all detected particles; 
Pz is the momentum along the beam axis and is the electron beam energy. 

Specific selection criteria applied in each decay channel are presented in the following sub- 
sections. A detailed description of the analysis can be found in [33]. 

5.1 Resonance Search 

The signature of the decay channel i/*— ^1/7 consists of an isolated electromagnetic cluster 
in events with missing transverse momentum. Background arises from CC DIS events with 
an isolated 7r° or a radiated photon. Events with substantial missing transverse momentum 

pmiss ^ 20 GeV are selected. In each event, a photon candidate with transverse momentum 
Py > 20 GeV in a polar angle range 5° < 9'^ < 120° is required. This polar angle range is 
restricted to 9^ < 60° in events with Pf"'^^ below 30 GeV, in order to reduce background from 
NC DIS. The photon is required to be isolated from jets by a distance R > 0.5 to any jet axis. 
In the central region (9'^ > 20°), photon candidates are selected only if no well measured track 
points to the electromagnetic cluster within a distance of closest approach (DCA) of 12 cm. For 
events with P™'^^ below 50 GeV, this condition is tightened by accepting only photon candi- 
dates having no track with a DCA to the cluster below 24 cm or within R < 0.5. The energy 
and polar angle of the photon are combined into one discriminant variable ^'^ = cos^ (9^/2). 
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Radiative CC DIS events are suppressed by requiring that > 45 GeV. For signal events, in 
most cases the final state contains a recoil jet, due to u* production via t-channel W boson 
exchange. Hence, in the final selection the presence of at least one jet with P;^^ > 5 GeV is 
also required. 

Seven events are selected in the data, compared to a SM expectation of 12.3 ± 3.0, which is 
dominated by CC DIS events. The invariant mass of the excited neutrino candidate is calculated 
from the four-vectors of the neutrino and the photon. The invariant mass distributions of the 
u* candidates and of the expected SM background are presented in figure [IJa). The resulting 
selection efficiency is 50% for M^. = 120 GeV, increasing to 55% for M^* = 260 GeV. From 
Monte Carlo studies, the total width of the reconstructed v* mass distribution is 11 GeV for a 
generated u* mass of 120 GeV, increasing to 41 GeV for a mass of 260 GeV. 

5.2 lyqq Resonance Search 

The signature of the u*^uZ^uqq decay channel consists of two jets with high transverse 
momentum in events with large P™"^*^. The SM background is dominated by multi-jet CC 
DIS events and contains a moderate contribution from photoproduction. Events with missing 
transverse momentum P^^^^ > 20 GeV are selected. In each event at least two jets are required 
in the polar angle range 5° < 9^'^^ < 130° with transverse momenta larger than 20 and 15 GeV, 
respectively. Additionally, the hadronic final state must exhibit a polar angle jh, as defined 
in [34], larger than 20°, in order to remove photoproduction events. Events with P™*^*^ < 30 GeV 
are selected if the ratio Vap/Vp of transverse energy flow anti-parallel and parallel to the hadronic 
final state [34] is above 0.1. This condition reduces the contribution of CC DIS processes. 
Photoproduction and NC DIS backgrounds typically have low volues of x^, the Bjorken scaling 
variable calculated from the hadronic system using the Jacquet-Blondel method [34,35], and 
are thus suppressed by requiring Xh > 0.04. Finally, to further reduce the background from CC 
DIS, a jet multiplicity greater than or equal to three is required for events with P™'^^ < 50 GeV. 
In each event, a Z candidate is reconstructed from the combination of the two jets with an 
invariant mass closest to the nominal Z boson mass. The reconstructed Z candidate is required 
to have an invariant mass above 60 GeV. 

After this selection, 89 events are found in the data compared to a SM expectation of 95 ± 21 
events. The invariant mass of the u* candidate is calculated from the neutrino and Z candidate 
four-vectors. For this calculation, the Z candidate four-vector is scaled such that its mass is set 
equal to the nominal Z boson mass. The invariant mass distributions of the u* candidates and 
of the expected SM background are presented in figure mb). The selection efficiency in this 
channel is 25% for M^* = 120 GeV, increasing to 55% for M^, = 260 GeV. From Monte Carlo 
studies, the total width of the reconstructed u* mass distribution is 31 GeV for a generated u* 
mass of 120 GeV, increasing to 41 GeV for a mass of 260 GeV. 

5.3 eqq Resonance Search 

The signature of the i'*^eW^eqq decay channel consists of one electron and two high Pt 
jets. Multi-jet NC DIS events constitute the main background contribution from SM processes. 
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Events are selected with an isolated electron in the LAr calorimeter in the polar angle range 
5° < 6''= < 90°. The electron variableQ = E^cos^ (6'72) is required to be above 23 GeV 
or the electron should have a transverse momentum greater than 25 GeV. These conditions 
remove a large part of the NC DIS contribution. In addition, the electron should be isolated 
from jets by a distance R > 0.5 to any jet axis. The events are required to have at least two jets 
in the polar angle range 5° < 9^'^^ < 130° with transverse momenta larger than 20 and 15 GeV, 
respectively. In each event, a W candidate is reconstructed from the combination of the two jets 
with an invariant mass closest to the nominal W boson mass. The reconstructed mass of the 
W candidate is required to be larger than 40 GeV. To further reduce the NC DIS background it 
is required that the polar angle of the jet with the highest Pt associated to the W candidate be 
less than 80° and that events with P^ < 65 GeV contain at least three jets with a Pt larger than 
5 GeV. 

After the selection, 220 events are observed compared to a SM expectation of 223 ± 47. 
The invariant mass of the u* candidate is calculated from the electron and W candidate four- 
vectors. For this calculation, the W candidate four- vector is scaled such that its mass is set 
equal to the nominal W boson mass. The invariant mass distributions of the u* candidates and 
of the expected SM background are presented in figure [TJc). The selection efficiency in this 
channel is 40% for M^* = 120 GeV, increasing to 65% for M^* = 260 GeV. From Monte Carlo 
studies, the total width of the reconstructed mass distribution is 15 GeV for a generated 
mass of 120 GeV, increasing to 38 GeV for a mass of 260 GeV. 

5.4 ez//i and eue Resonance Searches 

In the search for u*^eW—^eufj., events with P^^^^ > 15 GeV, one electron with P^ > 20 GeV 
and one muon with > 10 GeV are selected. The electron and the muon have to be detected 
in the polar angle ranges 5° < ^^"^ < 100° and 10° < 9^ < 160°, respectively. Furthermore, 
the electron and the muon must be isolated from jets by minimum distances of > 0.5 
and > 1, respectively. The contribution from NC DIS processes is reduced by requiring 
> 9 GeV. After this selection no data event remains, while 0.40 ± 0.05 SM background 
events are expected. The selection efficiency for v* with masses above 120 GeV is ~ 35%. 

The signatures of the i'*^eW^eiye and z/*— *z/Z— >z/ee channels are similar and consist 
of two high Pt electrons in events with large missing transverse momentum. Events with 
pmiss ^ 20 GeV are selected. In each event two isolated electromagnetic clusters are required, 
with a transverse momentum larger than 20 and 15 GeV, respectively. The highest Pt electron 
should be detected in the polar angle range 5° < 9^^ < 100° and the second electron in the 
range 5° < 9'^'^ < 120°. To reduce the background from Compton processes, a track has to be 
associated to each electromagnetic cluster in the central region (9^ > 35°). Events in which the 
invariant mass of the two electromagnetic clusters is within 10 GeV of the nominal Z boson 
mass are attributed to the Z— >ee decay channel. Events from the W^ve decay channel are 
selected by requiring the invariant mass of the u and one of the electromagnetic clusters to 
be compatible with the W boson mass within 20 GeV. In this channel, the variable defined 
from the highest Pt electron is required to be above 29 GeV. No data candidate is observed in 

' This variable is proportional to the four-momentum transfer squared for NC DIS. 
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either the ZorW decay channels compared to SM expectations of 0.19 ± 0.05 and 0.7 ± 0.1, 
respectively. In both channels, the selection efficiency for u* with masses above 120 GeV is 
~ 35%. 

5.5 Systematic Uncertainties 

The following experimental systematic uncertainties are considered: 

• The uncertainty on the electromagnetic energy scale varies between 1% and 3% depend- 
ing on the polar angle. The polar angle measurement uncertainty of electromagnetic 
clusters is 3 mrad. 

• The efficiency to identify photons is known with a precision of 10% for high Pt photons. 

• The scale uncertainty on the transverse momentum of high Pt muons amounts to 5%. 
The uncertainty on the reconstruction of the muon polar angle is 3 mrad. 

• The hadronic energy scale is known within 2%. The uncertainty on the jet polar angle 
determination is 10 mrad. 

• The uncertainty on the trigger efficiency is 3%. 

• The luminosity measurement has an uncertainty of 2.5%. 

The effect of the above systematic uncertainties are determined by varying the experimental 
quantities by ±1 standard deviation in the MC samples and propagating these variations through 
the whole analysis chain. 

Additional model systematic uncertainties are attributed to the SM background MC gen- 
erators described in section 3. An error of 20% on the simulation of NC DIS, CC DIS and 
photoproduction processes with at least two high Pt jets is considered to account for the uncer- 
tainty on higher order QCD corrections. An uncertainty of 30% on the simulation of radiative 
CC DIS events is considered to account for the lack of QED radiation from the quark line in 
the DJANGO generator. This uncertainty is estimated in the specific phase space of the anal- 
ysis by a comparison of the DJANGO result to the calculated cross section of the e'p^u^'yX 
process [36]. The error on the QED Compton cross section is estimated to be 10%. The errors 
attributed to multi-lepton and W production are 3% and 15%, respectively. The total error on 
the SM background prediction is determined by adding the effects of all model and experimental 
systematic uncertainties in quadrature. 

The theoretical uncertainty on the v* production cross section is dominated by the uncer- 
tainty on the scale at which the proton parton densities are evaluated. It is estimated by varying 
this scale from y^Q^/2 to 2\fQ^. The resulting uncertainty depends on the v* mass and is 10% 
at M^* = 100 GeV, increasing to 30% at M^* = 300 GeV. 
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6 Interpretation 



The event yields observed in all decay channels are in agreement with the corresponding SM 
expectations and are summarised in table[I] The SM predictions are dominated by NC DIS pro- 
cesses for the eqq resonance search and by CC DIS for the and uqq resonance searches. The 
distributions of the invariant mass of the data events are in agreement with those of the expected 
SM background as shown in figure [H No data event is observed in channels corresponding to 
leptonic decays of the W or Z bosons, in agreement with the low SM expectations. 

Since no evidence for the production of excited neutrinos is observed, upper limits on the u* 
production cross section and on the coupling //A are derived as a function of the mass of the 
excited neutrino. Limits are presented at the 95% confidence level (CL) and are obtained from 
the mass spectra using a modified frequentist approach which takes statistical and systematic 
uncertainties into account [37]. 

Upper limits on the product of the u* production cross section and decay branching ratio are 
shown in figure [21 The analysed decay channels of the W and Z gauge bosons are combined. 
The resulting limits on //A after combination of all decay channels are displayed as a func- 
tion of the v* mass in figure [3l for the conventional assumptions / = — /' and / = +/'. Limits 
are derived for u* masses up to 300 GeV. The total fraction of u* decays covered in this analysis 
is ~ 92% and ~ 84% in the two cases / = — /' and / = +/', respectively. In the case / = — /', 
the limit on //A is dominated at low mass by the z/*^z/7 channel, while the eW channel 
starts to contribute significantly for masses above 200 GeV. Under the assumption / = +/', the 
limit on //A is driven mainly by the u*^eW channel. These new results improve significantly 
the previously published limits by HI [2] and ZEUS [38]. For comparison, the most stringent 
limits obtained in e+e^ collisions at LEP for the two cases f = —f and / = +f', determined 
by L3 [39] and DELPHI [40], respectively, are also shown in figure [3l The HI measurement 
provides the most stringent constraints for masses larger than ~ 170 GeV. With the assump- 
tion //A = excited neutrinos with masses up to 213 GeV (196 GeV) are excluded for 
/ = -f (/ = +/'). 

Limits with less model dependence can be derived if arbitrary ratios /'// are considered. 
The dependence of the limits on this ratio for different u* masses is displayed in figure Ufa). 
Limits which are independent of /'// are derived for /'// G [—5; 5] by choosing in figureSfa) 
the point with the weakest limit for each mass hypothesis. The result is shown in figure |3b) 
and is found to be almost equal to the limit obtained under the assumption / = +/'. 



7 Conclusion 

Using the full e~p data sample collected by the HI experiment at HERA with an integrated 
luminosity of 184 pb^^ a search for the production of excited neutrinos is performed. The 
excited neutrino decay channels z/*^z/7, v*^uZ and v*^eW with subsequent hadronic or 
leptonic decays of the W and Z bosons are considered and no indication of a u* signal is found. 
New limits on the production cross section of excited neutrinos are obtained. Previous HERA 
results are improved by a factor three to four. Upper limits on the coupling //A as a function of 
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the excited neutrino mass are established for specific relations between the couplings (/ = +/' 
and / = — /') and independent of the ratio /'//. Assuming / = — /' and //A = 1/My*, 
excited neutrinos with a mass lower than 213 GeV are excluded at 95% confidence level. The 
results presented in this letter greatly extend the previously excluded domain and demonstrate 
the unique sensitivity of HERA to excited neutrinos with masses beyond the LEP reach. 
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Search for v* at HERA (e'p, 184 pb"^) 

Channel Data SM Signal Efficiency [%] 

v*^v-f 7 12.3 ± 3.0 50-55 

u*^eW^eqq 220 223 ± 47 40-65 

v*^eW^evn 0.40 ± 0.05 35 

v*^eW^eve 0.7 ± 0.1 45 

v*^vZ^vqq 89 95 ± 21 25-55 

v*^vZ^vee 0.19 ± 0.05 45 



Table 1: Observed and predicted event yields for the studied u* decay channels. The analysed 
data sample corresponds to an integrated luminosity of 184 pb~^. The error on the SM predic- 
tions includes model and experimental systematic errors added in quadrature. Typical selection 
efficiencies for v* masses ranging from 120 to 260 GeV are also indicated. 
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Search for v* at HERA (e p, 184 pb"") 
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Figure 1: Invariant mass distribution of the u* candidates for the z/*^z/7 (a), i'*^iyZ^iyqq 
(b) and u*^eW^eqq (c) searches. The points correspond to the observed data events and the 
histogram to the SM expectation after the final selections. The error bands on the SM prediction 
include model and experimental systematic errors added in quadrature. 
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Search for v* at HERA (e"p, 184 pb"*) 
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Figure 2: Upper limits at 95% CL on the product of the u* production cross section and decay 
branching ratio, ax BR, in the three decay channels as a function of the excited neutrino mass. 
The different decay channels of the W and Z gauge bosons are combined. Areas above the 
curves are excluded. 
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Figure 3: Exclusion limits at 95% CL on the coupling f/A as a function of the mass of the 
excited neutrino with the assumptions (a) / = — /' and (b) / = +/'. The excluded domain 
based on all HI e^p data is represented by the shaded area. Values of the couplings above the 
curves are excluded. The dashed line corresponds to the exclusion limit obtained at LEP by the 
L3 Collaboration [39] in (a) and by the DELPHI Collaboration [40] in (b). 
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Figure 4: (a) Exclusion limits at 95% CL on the coupling //A as a function of the ratio /'//. 
Each curve corresponds to a different u* mass. The circle indicates the weakest limit for each 
mass, (b) Exclusion limit at 95% CL on the coupling //A as a function of the mass of the 
excited neutrino (shaded area). This limit corresponds to the weakest limit on //A for /'// 
values in the interval [—5; +5]. The hatched area corresponds to the exclusion domain obtained 
by HI in a previous analysis [2]. The regions above the lines are excluded. 
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